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AFOSR-2278 

Abstract 

The  steady  flow  of  a  perfectly  conducting 
magnet ohydrodynamlc  fluid  past  a  thin  non-conducting 
airfoil  is  studied  in  the  model  in  which  the  fluid 
variables  obey  the  Lundquist  equations  linearized 
about  a  constant  unperturbed  flow.   "Hyperliptic " 
flows,  in  which  hyperbolic  and  elliptic  fields  are 
superimposed,  are  considered.   Results  of  Grad,  McCune 
and  Resler,  and  Sears  and  Resler  are  extended  to  the 
case  of  an  arbitrarily  inclined  unperturbed  field. 
The  general  solution  contains  four  line  singularities 
along  the  characteristics  through  the  ends  of  the  body, 
and  has  two  arbitrary  constants.   By  a  "generalized 
Kutta-Joukowski  condition"  these  constants  are  fixed 
so  that  two  of  the  line  singularities  disappear. 
Specifically,  it  is  required  that  the  solution  be 
locally  square  integrable.   Behavior  of  the  exponents 
of  the  singularities  is  investigated  by  numerical 
computation  and,  in  limiting  cases,  analytically.   The 
singular  parts  of  some  flows  are  investigated  numerically, 


-  2 


AFOSR-2278 

Table  of  Contents 

Page 

Abstract  2 

Section 

I.   Introduction  4 

II.   Formulation  of  the  Problem  6 

III.   The  Hyperliptic  Case:   General  Solution  13 

IV.   A  Generalized  Kutta-Joukowsky  Condition  21 

V.   Classification  of  the  Singularities  26 

VI.   Flow  Near  Singular  Edges  32 

Appendix  A  35 

Appendix  B  37 

Appendix  C  40 

Appendix  D  43 

Footnotes  48 

Figures  49 


3  - 


Magnet ohydrodynamic  Plow  Past  a  Thin  Airfoil 

I .   Introduction 

The  steady  flow  of  a  perfectly  conducting  magnetohydro- 
dynamic  fluid  past  a  thin  non-conducting  airfoil  has  been 
studies  previously. ^  '      The  standard  mathematical  model 
chosen,  and  the  one  used  here,  consists  of  the  Lundquist 
equations  linearized  about  constant  unperturbed  fluid  velo- 
city, density  and  magnetic  field,  together  with  the  appropriate 

(2) 
boundary  conditions.   '   Essentially  two  distinct  types  of 

flows  appear,  very  closely  corresponding  to  "subsonic"  and 

"supersonic"  flows  of  fluid  dynamics.   Since  the  regime  similar 

to  subsonic  flow  possesses  properties  common  to  both  subsonic 

and  supersonic  flows  the  situation  is  more  complex  than  either 

of  the  two  fluid  dynamic  cases.   This  case  is  the  more  inter- 

(2) 
esting,  and  with  the  methods  outlined  by  Gradv  ' ,    this  paper 

examines  the  "hyperliptic"  flows K    '    in  some  detail  and  extends 
the  results  to  cases  not  considered  previously.   In  much  of  the 
work  presented  all  vectors  are  considered  to  have  three  com- 
ponents, although  they  are  functions  of  two  space  variables 
only.   Some  sections  are  restricted  to  the  cases  where  vectors 
are  two  dimensional  and  where  the  unperturbed  flow  and  magnetic 
field  have  a  special  orientation. 

In  the  next  section  the  problem  is  formulated  and  Grad ' s 
results  are  given.  In  Section  3  a  precise  description  of  the 
hyperliptic  regime  is  given  and  the  formal  solution  to  the 
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equations  is  exhibited.   This  takes  the  form  of  a  representa- 
tion in  terms  of  Cauchy  integrals  over  the  body.   The  boundary 
conditions  on  the  body  then  yield  a  pair  of  coupled  integral 
equations  with  a  Cauchy-type  kernel,  whose  solution  is  known. 

The  solution,  however,  is  not  unique  and  in  Section  4 
an  uniqueness  condition  analogous  to  the  Kutta-Joukowsky 
condition  is  proposed.   '   Similar  to  subsonic  fluid  dynamics, 
functions  may  be  added  which  are  singular  at  the  ends  of  the 
body,  by  which  means  the  singularity  at  either  end  of  the  body 
may  be  suppressed.   The  dual  nature  of  the  hyperliptic  regime 
is  evident  in  the  fact  that  the  singularities  on  the  body  are 
also  propagated  into  the  fluid  along  the  real  characteristics, 
and  this  feature  is  used  to  determine  a  physically  acceptable 
uniqueness  condition.   The  condition  proposed  ensures  the 
square  integrability  of  the  flow  variables  over  any  region  of 
the  fluid.   The  line  singularities  remaining  in  the  solution 
result  from  a  breakdown  of  the  linear  analysis  used  and  the 
solutions  in  the  neighborhood  of  these  singularities  may  not 
provide  good  pointwise  approximations  to  the  correct  solution. 
On  the  other  hand,  they  may  give  an  indication  of  the  formation 
of  shocks. 

Since  the  exponents  of  the  singularities  are  complicated 
functions  of  the  flow  parameters,  numerical  values  were  obtained 
A  comprehensive  scheme  was  used  whereby  the  exponents  were 
found  over  the  whole  hyperliptic  region  for  a  set  of  flow  param- 
eters, and  the  regions  ever  which  the  trailing  and  leading 
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edge  singularities  are  to  be  removed  could  then  be  separated. 
This  is  described  in  Section  5.   In  the  strictly  two-dimen- 
sional problem  a  singularity  is  left  at  either  the  leading 
or  trailing  edge.   In  the  limit  of  vanishing  magnetic  field 
this  singularity  is  at  the  leading  edge  and  hence  the  ordinary 
Kutta-Joukowsky  condition  is  obtained.   In  the  general  problem 
a  singularity  is  left  at  each  edge.   The  exponent  of  the 
singularity  in  the  solution  of  the  restricted  problem  is  not 
obtained  as  a  uniform  limit  of  the  corresponding  exponents  in 
the  solution  of  the  general  problem^,  f  or  it  may  be  a  limit  of 
either  of  the  exponents  in  the  general  problem,  depending  on 
the  choice  of  U  .   The  flow  near  the  edges  of  the  body  was 
computed  for  several  cases  of  the  restricted  problem,  and  the 
results  are  given  in  Section  6. 


II.   Formulation  of  the  Problem 


The  steady  state  equations  for  a  perfectly  conducting 

plasma,  when  linearized  about  a  constant  flow  u  ,   magnetic 

field  B  ,   density  p  ,   sound  speed   a  ,   and  magnetic 
permeability  u.,   are 


(la)      (u  •  y)p  +  PnVu  =  0 


6  - 


(lb)      po(uQ-y)u  +  a2Q     p  -  i  (VXB)xB0  -  0 


(lc)       (uQ-  V)B  -  Vx(u  XBQ)  =  0 


(Id)       V-B  =  0   . 


Since  only  flows  depending  on  two  space  variables  are  con- 

rd   d 
^3x'  3y~ 


sidered,  the  gradient  operator  may  be  chosen  as  (*%—,    ^—  >    0) 


Further,  the  vectors  u  and  B  are  chosen  as 

o      o 


(2a)      uQ  =  uQ(l,  0,  0) 


(2b)       BQ  -  Eo(cos  6,    sin  6,    b)  , 


where  u   and  B  are  the  magnitudes  of  the  projections  of  u 
and  B   into  the  xy  plane.   As  written,  the  system  (l)  has 
eight  equations  but  only  seven  unknowns.   The  system  of  equa- 
tions actually  used  consists  of  (l)  with  the  x-component  of 
(lc)  replaced  by  the  condition  that  the  perturbed  flow  variables 
vanish  as  x  — >  +  co   for  y   fixed.   Under  this  condition  the 
omitted  equation  may  be  derived  easily  from  the  remaining 
equations . 

The  y-component  of  (lc)  may  be  integrated  immediately, 
and  with  the  same  assymption  on  behavior  at  infinity,  it 
follows  that 

(3)       B  cos  0  u   -  u  B   -  B  sin  0  u   =0. 
v  '  o       yoyo       x 
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Thus  u   may  be  eliminated  in  terms  of  u   and  B  .   With 
x    J  y       y 

the  introduction  of  non-dimensional  variables  defined  by 


(4a)  u  =  uQ(u,  v,  w) 


(4b)  B  =  B0(£,  n,  C) 


(4c)  P  =  pQa 


2     2 
and  on  the  definition  of  the  Alfven  speed,   AQ  =  Bo/u.po, 

the  basic  equations  become 


(5a)      Sin  9  ||  +  cos  9   *  +  sin  9  *  -  |l  =  0 


(5b)  a2sin   9  g  +   u2cos   9  ||  +    (A2sin29    -  u2)   §1 


.2     .    2„    d£  .  .2    .       Q    SC 

-   A    sin    0   -v2-  =    -    bA^sm    9   ^r 

o  ay  o  ox 

(50)  a2sin   9 (cos   9  ||  +   sin   9   *)   +  o2  ||  -  U2cos   9  & 

=    -   bA2sin   0(cos   0  |£-  +   sin   0  |£) 
o  v  ox  oy' 


(SO  H+&-° 


and 
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(6a)      u2  |^  -  A2 (cos  9   |£  +  sin  9   |£)  =  0 
v   ;        o  ox     ov       dx         dy' 

(6b)      sin  0(cos  0  |j  +  sin  0  |^)  -  sin  6   || 

i_  /     n  ^v        Q  dv    dn  \    „ 
-  b  cos  9  ^—  +   sm  0  -r—  -  ^-J  =  0 
v       dx         dy    dx' 

(and  u  sin  9  =   v  cos  9    -  r\ )  .      If  b  =  0,  the  variables  w  and 
^   disappear  from  (5)>  and  the  system  decouples  into  two 
independent  sets,  (5)  and  (6).   This  strictly  two-dimensional 
problem  will  be  considered  in  more  detail  in  some  sections 
below.   While  the  formulation  of  the  problem  is  incomplete 
without  boundary  conditions,  it  is  not  possible  to  assign  bound- 
ary values  without  a  study  of  the  solutions  of  the  systems  (5)? 
(6)  first. 

As  described  by  Grad,  the  system  (5)>  (6)  may  readily  be 
reduced  to  characteristic  form.   One  attempts  to  determine  a 
constant   A   and  a  linear  combination  of  the  variables   Z, 
Z  =  c-,o  +   CpV  +  c^w  +  oJc,   +  Cj-T)  +  Oft,    such  that  Z  satisfies 
the  equation 


^)  (S-X   fy)Z-0. 


Such  a  linear  combination  (called  a  Riemann  invariant)  exists 
if  and  only  if  A   satisfies  the  algebraic  equation 


-  9  - 


(8)  0   =    [(l+A2)(cos   0   +   A    sin   Q)2    -    (l+A2)(M2+  m2(l+b2)) 


o        o 


+   m2M2][(A    sin   9   +   cos    0)2    -   M2]     , 
o   o  o 


where  the  constants  m  and  M  are  defined  by  m  =  u  /a  , 

M  =  u  /A  .   For  the  real  values  of  A  satisfying  (8),  there 
o    o  o 

are  real  linear  combinations  Z  satisfying  (7)-   For  such  real 
A  (7)  states  that  Z  is  constant  on  the  characteristic  lines 
x  +  Ay  =  C .   The  complex  roots  of  (8)  generate  complex  linear 
combinations.   In  the  latter  case,  (7)  is  essentially  equiva- 
lent to  the  Cauchy-Riemann  equations  and  states  that  Z  is  an 
analytic  function  of  the  complex  variable  (x  +  Ay). 

The  second  factor  of  (8)  always  has  two  real  roots,  and 
the  nature  of  the  roots  of  the  first  factor  may  be  obtained 
by  the  simple  graphical  construction  described  by  Grad  (see 
Fig.  l).   Each  real  tangent  from  the  point  -u   to  the  charac- 
teristic cone  of  the  time  dependent  system  of  differential 
equations  generates  the  direction  of  the  characteristic  lines 
along  which  some  Z  is  constant,  and  a  positive  direction  along 
such  lines.   This  direction  is  from  -u  to  the  curve.   For  a 
given  characteristic,  the  positive  direction  is  called  "down- 
stream", because  it  is  the  direction  of  the  propagation  of  Z 
values  for  the  time  dependent  problem.   Correspondingly,  the 
negative  direction  is  called  "upstream".   The  terms,  "upstream" 
and  "downstream",  as  defined  here,  need  not  agree  with  the 
usual  definition  in  terms  of  direction  of  mass  flow.   However, 
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the  terms  are  relevant  in  that  it  is  at  the  upstream  ends  of 
characteristics  at  which  boundary  conditions  should  be  applied 
The  explicit  forms  for  the  various  Riemann  invariants  are 
given  by  Grad  and  are  exhibited  below.   It  is  convenient  to 
consider  a  vector  with  six  components, 


(9)        P  =  (^-,  v,  i,    T),  £,  w) 
mo 


and  represent  the  invariants  as 


(10)      z±  =   r±-P    i  =  1,  2,  3,  4,  5,    6. 


The  two  invariants  associated  with  the  quadratic  factor  of 
(8)  are  given  by 

(11)      r+  =  (l,-0,  +  ^75-  ,  +  -js-  ,  ^s-  ,  0)  + 

M      M    M 
o      o    o 

/tan20(l+A+)      \ 
[cos  0  t  MQ  -     -  ])l(0,-csc  0,0,0  +  ^-,  i) 

\ (1+a+  tan  6)  J  o 

where 


(12)       A+  =  (±  MQ  -  cos  0)cse  0. 


The  four  other  invariants  may  be  expressed  conveniently  in 
terms  of  the  parameters  defined  by: 


1    In 
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(13)  tan  j6  =  A 

(14)  6  =  9    -  j6    , 


so   that 


(15)  cos25.     -   cos2$.[M2   +(l+b2)m2]    +   mV    cos  j6±   =   0, 

i  =   1,   2,    3,    4 
and 


(16)       r.  =  (M  sin  j6. cos  j6 .     -    sin  6   cos  5^  cos  ^(l-lY^cos  j6± )  , 


-  cos  J#.sin5..   sin  (#.cos  f^.sin  6. 


b  cos  $.sin  j6.     ,      b  cos  j^.sin  0) 


Consistent  with  the  linearization,  the  boundary  conditions 
will  be  imposed  on  a  slit,  for  convenience  taken  as  y  =  0, 
|  x|  <  1.   The  limit  of  a  quantity  taken  as  a  point  approaches 
the  top  or  bottom  of  the  airfoil  is  denoted  by  the  subscript 
t   or  b   The  condition  that  the  velocity  normal  to  the  body 
vanishes  at  the  body,  when  linearized,  gives  v   and  v,  as  known 
functions  of  x.   It  is  also  assumed  that  the  magnetic  field 
is  continuous  across  the  body  and  further,  from  the  condition 
V  x  B  =  0  inside  the  nonconducting  body,  £.  =  £,  =  constant. 
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The  remaining  conditions  depend  on  the  nature  of  the  solutions 
considered.   All  "eliptic"  Riemann  invariants  are  assumed  to 
vanish  at  infinity  and  no  further  restriction  is  imposed.   For 
the  hyperbolic  invariants  the  situation  depends  on  the  orien- 
tation of  the  characteristics.   If  the  characteristic  does  not 
pass  through  the  body,  then  the  invariant  vanishes  on  that 
line.   For  characteristics  passing  through  the  body,  the 
invariant  vanishes  on  the  half-line  which  has  infinity  "up- 
stream".  On  the  other  half-line  no  condition  is  imposed. 


Ill .   The  Hyperliptic  Case:   General  Solution 

In  the  region  of  interest  in  this  paper  the  character- 
istic equation  (8)  has  one  pair  of  complex  conjugate  roots, 
denoted  by  A?  =  s  +  it   and   A  -  A  .   The  Riemann  invariants 
Z.,   i  =  1,  4,  5>  6   corresponding  to  the  real  characteristic 
roots,  are  zero  on  half  lines  leaving  the  body  and  with 
infinity  upstream.   Thus,  by  convention  the  boundary  condi- 
tions may  be  taken  as 


(17a)      (Z±)t  =0      i  =  1,  5        on  |x|  <  1 


(17b)     (Z1)b  =0      i  =  4,  6  . 
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For  the  complex  roots  of  (8),  if  Z  is  defined  by 

(18)       z  -  x  +  sy  +  ity  =  X  +  iY  , 

then  Zp  is  an  analytic  function  of  z,    which,  by  the  boundary 
conditions,  vanishes  at  infinity.   Under  the  transformation 
implied  by  (l8),  the  slit  is  unchanged 

The  function  Zp  may  be  expressed  in  terms  of  its  real 
and  imaginary  parts  as 


Z2  =  Wl  +  1W2 


Since  Z   is  an  analytic  function  of  z  =  X  +  iY,  it  may  be 
represented  as  a  Cauchy  integral  in  terms  of  its  boundary 
values.   Under  the  assumption  that  Zp  tends  to  zero  at  infinity 
sufficiently  rapidly,  this  representation  has  the  form 

(19)     Zp  =  k  r  p(p)  +  iq(p)  dp  . 

v  J '  2    n  J      p-z       v 

-1 

Thus  Zp  must  go  to  zero  at  infinity  at  least  as  fast  as  z 
The  functions  F  and  G  in  (19)  are  defined  by 

(20a)     F(X)  =  \   {  (W1)t  -  (W1)bj 

(20b)      G(X)  =  \\   (W9),  -  (W9), 


2   \"2't    v "2'b 
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The  boundary  values  of  w\  ,  W  are  determined  from  (19)  as 

1 
(21a)     (W1(X))t  =  P(X)  +  i  f  f-^dp 

J-l 


(21b)      (W1(X))b  =  -P(X)  +  i 


'  Gl£] 


7T 


dp 


(21c)      (W2(X))t  =  G(X)  -  I  [  |i£) 


J 


1 
1 

dp 
-1 


1 

(21d)      (W2(X))b  =  -G(X)  -  \   [   ^-dp 

-1 


where  the  integrals  are  defined  as  Cauchy  principal  values. 
Clearly  (20)  and  either  equation  for  W-.  (X)  and  for  Wp(X)  from 

(21)  is  equivalent  to  the  system  (21).   Thus  the  fundamental 
equations  become  (17)?  (20)  and  two  appropriate  equations 
from  (21). 

The  system  of  equations  (17)  and  (20)  may  be  considered 

a  system  of  six  linear-  equations  in  the  eight  unknowns  o,,    a,  , 

£>  r\  y    w,  ,    w,  ,  F,  and  G.   The  inhomcgeneous  terms  in  the 

system  are  linear  functions  of  v. ,  v,  and  t    ,    functions  given 
J  t   b      o  D 

by  the  boundary  conditions.   Thus  o,  ,    o,,    £, ,  i],    w,  ,  w,  may  be 

expressed  as  a  linear  combination  of  F,  G,  and  the  known 

functions  v, ,  v,  and  £    .      The  form  of  the  two  integral  equa- 
t    b      bo  e      m 

tions  taken  from  (21)  becomes 

(22)  AH  +  IB  J  ^"dp  =  C 

-1 
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where  A  and  B  are  two-by-two  matrices  with  constant  coeffi- 
cients, and  C  is  a  two -component  vector  which  is  a  known 
function  of  X  (a  function  of  v,  ,  v,  and  £  )  and  H  represents 
the  two-component  vector  (F,G).   The  matrix  A   is  assumed  to 
be  nonsingular,  and  pre-multiplication  of  (22)  by  A   gives 


(23) 


1   i 

H  +  -  B 

IT 


-1 


|Mdp=c. 


There  exists  a  nonsingular  matrix  Q,  such  that  B  =  QB'Q 
is  in  sub -diagonal  form 


-1 


(24) 


where  a  £   0  only  if  p,  =  v.   Under  the  transformation  H  =  QH 
equation  (23)  becomes 


(25) 


3  +  ^1    f^^p 

J-l 


QC   , 


which  represents  the  pair  of  ihhomogeneous  singular  integral 
equations  with  a  Cauchy-type  kernel 


(26a) 


(26b) 


Hn    +   ^ 

1  IT 


J 


r1  Mp) 
-jhr  dp  =  (QC  }i 


H2    +   7 


-1 
1 

"-1 


h2(p)  ,         a  r    h  (P) 

dp  =  (qC) 


p-X 


2  7T  p-X 

-1 


dp 
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2    ?  /x 

Provided   |i  ,  v  ^  -1 ,  equation  (26a)  is  solved  for  H-,  by 

standard  techniques  (see,  for  instance,  Mikhlin^-^  )  and  the 

solution  is  substituted  in  the  right  hand  side  of  (26b)  which 

is  then  solved  for  H0 . 

The  solution  of  an  integral  equation  of  the  type 

1 


(27) 


( 


•(X)  +  % 


MP) 

p-X 


dp  =  f(x) 


is  given  by 


(28)         fi(x; 


1+5 


f(X) 


r 


1   5    /-,,-.,  \m-l  /..  „v-m  I   ,.,   xl-m,,   \mw  w   -u-\~l  j 
-  p"  (1+X)    (1_X)    I   ^1+P)    ^1-P)  f(p)(p-X)    dp 


-1 


+  Kd+D^d-x)1"  , 


where   m   is  determined  from 


(29) 


2  Trim    1+15 
~  1-15 


0  <  Re(m)  <  1 


It  is  noted  that  the  solution  (28)  contains  the  arbitrary 
constant   K.   If   f(x)   is  well-behaved,  e.g.,  Lipshitz  con- 
tinuous, then  it  may  be  shown  easily  that 


(30a) 


j6(x) 
x  — >  1 

x  <  1 


1 

,1-m 


(C+K) 
( .1  -x ) 


0((l-x) 


1-m, 
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(30b)  *S(x)      =1--     (C'fflm+0((l+x)m)    , 

x  —>-l        2         (1+x) 
x  >    -1 


where 

1 

(51a)  C  = ^-   I    dp(f(p)    -  f(l))     ^g-i 

tt(1+5    )  J  \      H/ 

-1 

1 

(31b)        c  =  -77^7  I  dp(f(p)  -f(-1))    fep 

"-1 


Hence   K  may  be  chosen  so  that  j6(x)    Is  regular  (and  even 
vanishes)  at  the  leading  edge  or  at  the  trailing  edge.   However, 
j6(x)   will  be  singular  at  one  edge  or  the  other  unless   f   is 
specially  chosen.   Such  a  choice  of   f   implies  a  special 
shaping  of  the  body,  as  is  considered  briefly  below  (p.  24). 

With  the  solutions  to  the  integral  equations  of  the  form 
(26),  it  remains  to  describe  the  solution  of  the  original  prob- 
lem.  H-,  (x)  and  H?(x),  linear  combinations  of  F(x)  and  G(x), 
satisfied  (26);  thus  H,(x)  and  Hp(x)  may  be  written  in  the  form 
(28)  with  constants  m-,  ,  K,  and  m? ,  Kp .   Further,  near  the  ends 
of  the  body,  H, (x)  and  Hp(x)  have  the  behavior  given  in  (30) 
with  the  appropriate   m  and   K.   It  follows  that  F(x)  and 
G(x)  are  each  the  sum  of  two  expressions  of  the  form  (28)  with 
constants  rru  ,  K,  and  m  ,  Kp.   Since  Zp  ( z )  is  given  by  (19)  >    it 
is  easily  seen  that  for   z   near  the  leading  or  trailing  edge, 
i.e.,  y  — ^  0,  x— t-+  1, 


18 


(32a)         Z  (z)   =    a  m  +    %p  +  0(l) 

(1-z)  l    (1-z)  2 


z— >>  +1 


(32b)         Z2(z)   =- 7I3mT"+- ^2  +  0(1). 

z_^  _!    (l+z)    -1    (1+z)    ^ 


Again,  with  judicious  choice  of  the  two  constants   K-,   and 
K0 ,  one  constant  of  each  pair  of  constants  (a, a')  and  ((3,p') 
may  be  set  to  zero. 

Now  that  the  elliptic  part  of  the  solution  has  been  com- 
pleted, the  hyperbolic  part  must  also  be  given.   Since   F(x) 
and  G(x)  are  known,  equations  (17)  and  (20)  constitute  six 
equations  for  the  six  unknowns  a ,  ,  a    ,  £, ,  t\,   w ,  ,  and  w  . 
Thus  these  variables  are  given  as  a  linear  combination  of 
P(x),  G(x)  and  the  known  functions  v, (x),  vK(x),  and  £  .   The 

Xf  0  o 

Riemann  invariants  previously  unspecified  on  the  body  have 
the  form 


(33a)         (zi(x))b  =  aip(x)  +  PiG(x)  +  V'x)' 

Ixl  <  l,  i=i,5; 


(33b)         (Z1(x))t  =  a±P(x)  +  P±G(x)  +  R±(x), 

Ixl  <  1,   i  =  4,  6  , 


where  R. (x)  are  well-behaved  functions  of  x.   Since  the  i 
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Riemann  invariant  is  constant  on  the  line   x  +  yA. ,  the  non- 
zero values  of  the  Riemann  invariants  are  given  everywhere  by 
(33)  on  the  replacement   x  — >•  x  +  yA.  .   Hence  the  singularities 
present  in  F(x)  and  G(x)  at  the  tips  of  the  body  propagate  along 
characteristics.   In  general,  the  real  Riemann  invariants 
behave  like 


(34a)     Z  (x,y)   =- 1 -*=-+- T"   m9  +  0(1) 

■\        . -,    (l-x-yA.   -1    (l-x-yA.)  ^ 

x+yAj_  <  1 


(34b)     z±(x,7)   =  zLT35-h-  ^^-T^  +  0(D  • 

x+yA±— ^-1    (l+x+yA1)        (l+x+yA±) 
x+yAi  >  -1 


The  choice  of  K-,   and  Kp  which  removed  a  singularity  of  a 
complex  invariant  eliminates  the  corresponding  singularities 
in  the  real  invariants.   The  value  of  any  Riemann  invariant  on 
its  characteristic  is  zero  unless  the  characteristic  ends  on 
the  body,  and  then  its  value  is  given  by  (17)  or  (33).   Since 
each  physical  variable  is  a  linear  combination  of  all  the 
Riemann  invariants,  each  variable  has  the  composite  behavior 
implied  by  (32)  and  (3*0 .   Hence,  as  one  approaches  a  char- 
acteristic leaving  a  tip  the  behavior  anywhere  in  the  fluid 
is  given  by  (3^).   As  one  approaches  a  tip  itself  the  behavior 
is  the  sum  of  terms  from  (32)  and  (3^). 

An  expression  of  the  form  l/|l-x-yA|    is  square - 
integrable  on  line  segments  or  over  areas  only  if  m  <  l/2. 
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Hence,  the  singularities  in  (3^a)  and  (3^b)  from  a  given 
exponent  are  not  both  square-integrable.   The  singularities  in 
the  complex  invariant,  see  (32),  are  square-integrable  over 
areas  provided  0  <  m  <  1.   The  generalized  Kutta-Joukowsky 
condition  of  the  next  section  is  based  on  these  observations. 

In  the  case  b  =  0,  the  variables   £   and  w   separate 
out  and  the  Riemann  invariants  Z,-  and  Zr     may  be  ignored 
or  considered  separately.   Now  it  is  possible  to  obtain  a 
linear  relation  between  F(x),  G(x)   and  known  functions. 
Hence  there  is  only  one  integral  equation  and  one  exponent 
m  and  constant   K.   Otherwise  there  are  no  changes  from  the 
above  description. 

IV.   A  Generalized  Kutta-Joukowsky  Condition 

The  present  section  deals  with  the  choice  of  the  constants 
K-,  ,  K„  necessary  for  a  unique  solution.   In  ordinary  fluid 
dynamics  a  similar  non-uniqueness  occurs  in  the  subsonic  solu- 
tion:  in  this  case  a  single  integral  equation  is  solved  m  =  p-, 
and  the  solution  contains  one  arbitrary  constant.   The  Kutta 
condition  requires  smooth  flow  at  a  sharp  trailing  edge  and  is 
equivalent  to  the  addition  of  circulation  to  the  flow  sufficient 
to  cause  a  confluence  of  the  rear  stagnation  point  on  the  body 
with  the  trailing  edge.   The  singularity  at  the  leading  edge 
remains  in  the  linearized  solution  in  general  unless  the  body 
is  shaped  to  remove  it. 

A  generalization  of  the  Kutta-Joukowsky  condition  to 
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determine  a  unique  solution  in  the  magneto -flu id  dynamic  case 
is  proposed  here.   The  singularities  in  this  case  have  expo- 
nents m  and  1-m,  (0  <  m  <  1),  at  the  trailing  and  leading 
edges,  respectively,  where   m  may  have  two  values  in  the 
general  two-dimensional  problem.   The  proposed  condition  is 
to  require  the  square  integrability  of  these  variables  over 
any  finite  region  in  the  flow,  in  particular  in  the  neighbor- 
hood of  the  line  singularities.   This  condition  uniquely 
specifies  the  two  constants  K-,  and  K? .   In  an  exact  solution 
such  square  integrals  must  exist  as  they  represent  the  energy 
in  a  region  or  the  forces  acting  there.   If  the  proposed 
condition  were  not  valid,  then  the  linearized  solution  would 
be  a  bad  approximation  to  the  exact  solution  near  the  line 
singularities  both  pointwise  and  in  the  square  integral  sense. 
Thus,  non-linear  or  dissipative  phenomena  would  affect  the 
gross  characteristics  of  the  flow  and  the  linearized  solution 
would  have  almost  no  meaning.   Hence,  if  the  linearized  solution 
is  to  be  taken  as  a  valid  approximation  (except  in  certain 
regions)  then  it  must  be  used  in  conjunction  with  the  proposed 
condition.   The  weaker  point  and  line  singularities  remaining 
in  the  flow  are  not  expected  to  be  good  pointwise  approximations 
to  the  correct  flow  there,  and  a  non-linear,  dissipative  analysis 
may  have  to  be  done  to  improve  the  solution  in  these  regions. 
Nonetheless,  one  might  hope  that  the  over-all  characteristics 
of  the  linearized  solution  are  reasonably  correct,  as  in  ordinary 
fluid  dynamics  where  integral  properties  such  as  lift  on  the  air- 
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foil  are  predicted  fairly  accurately  by  the  linearized  theory. 

However,  experimental  results  will  be  the  final  arbiter. 

For  the  restricted  problem  a  curve  in  the  hyperliptic 
region  was  found  which  separates  regions  m  >  p-  from  m  <  p-.  On 

the  two  sides  of  this  critical  curve,  m  has  the  values  1  and  0. 
It  is  instructive  to  follow  the  behavior  of  the  singularities 
as  the  critical  curve  is  crossed,  say,  from  a  region  m  <  p- 
to  a  region  m  >  p-  (e.g.,  by  varying  u  ).   In  the  region  m  <  p- 
the  leading  edge  singularity  is  removed.   As  the  critical  curve 
is  approached  the  trailing  edge  singularity  tends  to  zero  and 
at  the  critical  curve  there  is  no  singularity  at  either  edge. 
On  entering  the  region  m  >  p-,   m  is  close  to  unity,  the 
trailing  edge  singularity  is  removed  and  the  leading  edge 
singularity  is  small.   The  total  effect  of  crossing  the  critical 
curve  is  that  the  singularity  disappears  from  one  edge  and 
reappears  at  the  other.   It  is  noted  that  this  is  a  continuous 
process  and  Is  dependent  on  the  Kutta  condition  used.   If,  say, 
a  condition  that  the  trailing  edge  singularity  be  removed  were 
used  instead,  then  there  would  be  a  discontinuous  disappearance 
of  the  singularity  at  the  leading  edge  across  the  critical  curve. 

A  Kutta  condition  having  similar  effect  to  the  one  des- 

(6) 
cribed  above  has  been  obtained  by  Stewartson^  '    for  the 

restricted  problem  in  the  incompressible  limit.   The  removal 

of  one  of  the  edge  singularities  in  either  region  m  >  p-  or  m  <  p- 

is  related  to  its  removal  on  the  relevant  side  of  the  critical 

curve.   There, the  singularity  at  one  of  the  edges  has  the  limit 
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unity  for  its  exponent  and.  when  m  =  0   (or  l)  there  are 
infinite  forces  on  the  body  in  the  neighborhood  of  the  leading 
(or  trailing)  edge.   Since  such  a  solution  is  unacceptable, 
Stewart son,  by  a  continuity  argument,  removes  that  singularity 
whose  exponent  tends  to  one  as  the  critical  curve  is  approached 
from  the  relevant  side.   Thus  Stewartson's  condition  produces 
the  same  result  as  the  one  proposed  here,  but  it  is  imposed 
for  different  reasons.   The  argument  here  is  based  on  the 
unacceptability  of  the  solution  without  the  proposed  Kutta- 
Joukowsky  condition  without  recourse  to  any  continuity  argu- 
ments.  The  distinction  is  clear  in  the  general  two-dimensional 
case  where  Stewartson's  criterion  is  inapplicable  as  computa- 
tions have  shown  that  the  two  exponents  m,  and  m2  never  become 
zero  or  one . 

The  singularities  considered  previously  are  those  arising 
in  the  solution  to  the  homogeneous  integral  equation.   The 
shape  of  the  body  near  the  leading  and  trailing  edges  also 
determines  the  singular  parts  of  the  solution.   For  convenience, 
the  edge  of  the  body  is  now  taken  at  the  point  x  =  0  and  the 
body  oriented  to  lie  in  X  >  0.   In  a  neighborhood  of  x  =  0 
the  body  is  assumed  to  be  represented  on  top  and  bottom  by 
the  functions  y  =  h^x),  h2(x)  of  the  form  hi(x)  =x(c±1  +  xR±1  (x) ) 
+  xr(c2  +  xR.p(x))   where   0  <  r  <  1   and  R.  .   are  regular 
functions  of  x.   If   c.p  4-   0  the  term  c.pXr  dominates  at 
an  edge  and  represents  a  blunt  body. 

The  inhomogeneous  term  in  the  integral  equation,  depending 
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linearly  on  the  slope  of  the  body,  has  the  form 


f(x)  =  d1(l  +  xR1(x))  +  xr  Xd2(l  +  xR2(x)) 


where  R. (x)   is  regular.   The  constants  d.  are  scalars  or 
matrices  according  to  whether  a  single  equation  or  a  system 
is  considered.   With  ex   representing  the  singular  part  of 
the  non-trivial  solution  of  the  homogeneous  equation,  the 
solution  of  the  inhomogeneous  equation  near  the  given  edge 
is  of  the  form 


e-,/x  +  eyx    +  R(x)  , 


■y  >    -2< 


where  R  is  regular  and  e-.  is  determined  by  the  generalized 
uniqueness  condition.   If   m  <  p-,   e-,  is  determined  by 
application  of  the  condition  at  the  other  edge  of  the  body. 
If   m  >  _-,   e-.  is  put  equal  to  zero.   In  either  case,  the 
solution  contains  a  singularity  of  order  (l-r).   For  the 
above  arguments  in  favor  of  the  proposed  uniqueness  condition 
to  be  valid,  it  must  be  assumed  that   r  >  p-   (the  case 
l-r  =  m  >  p-  is  dismissed  as  too  coincidental),  and  hence 
some  blunt  bodies  are  admissible.   If  c,?  =  0  (i  =  1,  2)  so 
that  the  body  is  not  blunt,  the  term  in  l/x    is  absent. 
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V.   Classification  of  the  Singularities, 

The  exponent   m  of  the  singular  terms  involved  in  the 
solution  of  the  integral  equations  (26)  are  complicated  func- 
tions of  the  unperturbed  velocity  and  magnetic  fields  u  ,  B 
and  the  Alfven  and  sound  speeds  A  ,  a  .   Numerical  calculations 
were  performed  to  investigate  the  variation  of  m  with  these 
parameters . 

In  the  restricted  two-dimensional  problem,  m  has  a  single 
value  which  may  be  obtained  directly  from  the  coefficient  5  of 
the  integral  equation  corresponding  to  (27)  for  F  or  G  (recall 
that  P  and  G  are  related).  The  formulae  necessary  for  this  cal- 
culation are  set  out  in  Appendix  A.  Computations  were  made  to 
find  m  in  the  hyperliptic  region  (that  is,  over  the  whole  range 

of  values  of  u   and  6)    for  a  range  of  values  of  the  ratio 

2   2  1     -1 

a  =  a  /A  .   When  the  coefficient  5  is  positive,  m  =  —  tan   6, 

and  0  <  m  <  -p.   When  5  is  negative,  m  -   1  +  —  tan   6,  and 
■p  <  m  <  1.   Hence,  where  5  changes  sign,  m  jumps  from  zero  to 
unity;  this  curve  separates  regions  of  m  <  p-  from  m  >  p- 
and  is  called  the  critical  curve. 

Figures  2  to  7  indicate  the  variation  of  the  critical  curve 
in  the  restricted  problem  over  the  range  of  values  of  a,  and 
various  contour  lines  of  m  are  also  shown.   The  figures  represent 
one  quadrant  of  the  hyperliptic  region  and  the  complete  diagram 
can  be  obtained  by  reflection  in  the  axes.   The  unperturbed  mag- 
netic field,  B  ,  is  taken  along  the  horizontal  axis  (as  in  Fig. 
1).  The  velocity  vector,  u  ,  extends  from  the  origin, and  the 
size  and  inclination  of  this  vector  with  respect  to  the  other 
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parameters  then  determine   m  and  on  which  side  of  the  critical 
curve  it  lies.   Flows  with  u   and  B  perpendicular  always  lie  in 
the  region  m  >  „-.   For  large   a   the  critical  curve  is  a  con- 
tinuous curve  from  the  origin  to  the  cusp  on  the  axis  of  the 
inner  wave  front.   The  region  m  <  p-  increases  in  area  as  a 
decreases  and  splits  into  two  branches  as  it  overlaps  with  the 

outer  limit  of  the  hyperliptic  region  (Figure  h) .      The  branch  to 

2    2 
the  inner  wave  front  is  lost  as  the  limit  a  =  A   is  passed  and 

o    o    ^ 

for  a  <  1  the  critical  curve  extends  from  the  origin  to  the 

/"" 2    2"  

a  +A  )  on  the  perpendicular  to  B  .  Figure  2,  corres- 
ponding to  the  incompressible  limit  (A  fixed,  a  — >  <») ,  shows 
Stewart  son's  result  that  the  critical  curve  is  given  by 

(35)  M2  =  cos  20  . 

Alternatively  it  represents  the  limit  of  vanishing  magnetic  field 

(A  — ^  0,  a   fixed)  in  which  case  the  m  <  p-  region  is  very  small 

and  the  Kutta-Joukowsky  condition  of  ordinary  fluid  dynamics 

is  recovered  as  a  limit  of  the  generalized  uniqueness  condition. 

The  results  just  obtained  may  be  checked  with  information 

inferred  from  certain  expansions  of  5.   For  small  values  of  u 

^  o 

the  coefficient  5  of  the  integral  equation  in  the  restricted 
case  has  the  behavior 

(36)  5^am  cos  29 /[sin  9  /l+a  }   . 

Thus  the  critical  curve  has  slope  -jr     at  the  origin.   In  the 
limit  9  -^  0,      u   arbitrary  but  fixed, 
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(37)  5-  2amo(l-m2)(_L_  -  m2)/jec/l+o7  [a2m2+(l+a) (l-m2 )2 ] j 


where 


,  l/2 
(38)  C   =   /(l-m2)(l-M2)/(l-   2    -  M2)i 


Hence   5   changes  sign  twice,  once  at  u  =  a  and  once  at 

2    2* 
a  +  A   ,  the  latter  point  being  the  position  of 
00  ^ 

the  cusp  of  the  inner  wave  front  on  the  B   axis.   For  a  >  1 

(and  hence  9   =   0,  u   =  a   lies  on  the  outer  wave  front),  6  has 
v  00 

opposite  signs  in  the  hyperliptic  region  near  9   =   0  on  the  dif- 
ferent sides  of  the  cusp  and   6  vanishes  at  the  cusp.   As  was 
observed  in  the  figures,  the  critical  curve  leaves  the  origin 
and  reenters  the  cusp.   For  a  <  1  (and  hence  0  =  0,  u   =  a 
lies  on  the  inner  wave  front),  in  the  hyperliptic  region  near 
9   =   0,  6  has  only  one  sign.   Now  the  critical  curve  leaves  the 
origin  and  terminates  at  8  =  ■£,   u   =  /a  +  AQ  .   Sears  and  Grad 
have  shown  that  when  u   and  B  are  parallel  the  system  behaves 
as  an  ordinary  fluid,  while  from  (36)  as  9   — 7  0,   m  — ^  +  75", 
the  fluid  dynamic  result.   Without  much  more  care  little  more 
can  be  inferred  from  this  treatment  of  the  parallel  flow  case. 
Numerical  calculations  to  determine  the  exponents  m.  ,  iiu 
in  the  general  two-dimensional  case  were  also  performed.   The 
equations  leading  to  the  coupled  integral  equation  (20)  are 
given  in  Appendix  B.   The  numerical  results  indicate 
0<m1<p-<m2<l,   where   m.  —?  75-  as  9  — ^  0. 
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According  to  the  proposed  uniqueness  condition,  there  is  a 
singularity  of  exponent   rru   at  the  leading  edge  and  a  singularity 
of  exponent   1-m   at  the  trailing  edge.   Thus  the  leading  edge 
singularity  is  stronger  or  weaker  than  the  trailing  edge  one  as 
m.,  +  m   is  greater  or  less  than   1  .   In  figures  9-15  the 
curves   m.,  +  m  =  1   are  exhibited  for  various  values  of  the 
parameters.   The  unit  of  speed  is  the  Alfven  speed  A   .   For  b 
large  the  two  branches  of  the  critical  curve  do  not  fit  conveniently 
on  one  scale,  hence  part  of  the  figure  is  shown  in  a  reduced 
scale  on  the  left.   In  such  cases  the  second  branch  is  supressed 
in  the  right-hand  figure.   The  inner  wave  front  either  is  drawn 
or  if  it  is  too  small  its  position  Is  indicated  by  a  square. 
In  each  figure  one  region  is  marked  for  which  iru  -t-  mp   is 
greater  or  less  than  one.   In  all  cases,  as  the  curves  are  crossed, 
the  inequality  changes. 

For  small  values  of   b   one  exponent  should  be  close  to  the 
one  computed  when  b  =  0  and  the  other  should  be  near  zero  or  one. 
Hence  the  curve  on  which  m  =  0   ( or  1 )  for  b  =  0   should  ap- 
proximate the  curve  m,  +  mp  =  1   for   b  =f=  0  but  small.   Thus 
Figures  9  and  13  are  similar  to  Figures  6  and  3  respectively. 
However  in  each  case  there  is  one  additional  curve  when  b  f  0  , 
but  this  curve  vanishes  as   b   goes   to  zero.   Hence  the  distinc- 
tions present  between  a  >  1   and  a  <  1   when  b  =  0  also 
occur  for  small  values  of   b  =(=  0  .   As  shown  in  the  figures  10, 
11  and  12,  as   b   increases  the  nature  of  the  curves  changes  for 
a  <  1  ,  so  that  the  distinctions  between  a  <  1   and  a  >  1 


-  29  - 


disappear. 

Because  of  various  difficulties  in  the  numerical  work  certain 

features  of  figures  9-15  cannot  be  readily  verified.   First,  it 

is  not  definitely  known  that  the  curves  enter  the  origin,  as 

shown  .   Secondly,  it  is  not  clear  that  the  curves  shown  to 

intersect  the  outer  surface  of  the  ray  cone  do  not  continue  very 

close  to  the  cone  and  end  at  some  other  place.   In  both  cases  the 

numerical  work  strongly  indicates  the  conclusions  implied  by  the 

figures . 

For  large   b  ,  and  m   and  M   fixed,  the  coefficients  \±,v 

o        o 

of  the  integral  equation  (26)  are  given  by 

M  (cos  29  -  M2) 
(39a)  M.--2 J2_  +  o(I) 

sin  e    (1  +  ]Yr)      " 

o 

-m  M2  sin  0  . 

(39b)  v-  — ^ 5-  +  0(±o)   . 


b  (cos  20  -  ivp)      b 


2 


An  alternate  description  of  the  critical  curve  would  be  a  curve 
across  which  \i  +   v   changes  sign.   From  (39a)  and  (39b)   one 
branch  of  the  curve  for  b   large  is  seen  to  be 


2 


(40)  M  =  cos  20 


o 


O  Q  O 

The  other  branch,  found  by  setting  m   and  M  =   0(b  )   ,  was 
computed  to  be 


(4la)  M2  +  m2  sin20  =  b2 

v    '  00 
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2      2 
or,  with  M  =  am  , 
o     o 


(4lb)  M2 


o        .  2n 

a  +   sin  9 

2         2       2 
For  m   and  Mo  =  0(b  ),  the  coefficients   u.  and   v   of 


the  integral  equation  (26)  were  found  to  be 


m   sin   6        m   sin   0 
(42a)  n~ 


2 
am 

<42b)  v~   m   sln"e   =  '  m   sin   6 

o  o 

2   2      2   2 
where   M  /b   =  am  /b  <  1   is  necessary  if  the  velocity  u   i 

to  remain  in  the  hyperliptic  region  for  large   b.   Without 

further  terms  in  the  expansion  (4l)  it  is  not  possible  to 

decide  whether  this  branch  of  the  critical  curve  intersects 

the  characteristic  cone  for  large   b. 


s 
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VI.   Flow  near  singular  edges. 

In  order  to  understand  the  nature  of  the  singularities 
present  in  this  problem   the  flow  quantities  in  the  neighborhood 
of  the  leading  and  trailing  edges  have  been  computed  numerically. 
The  flows  are  also  exhibited  for  singularities  rejected  by  the 
proposed  uniqueness  condition.   The  flows  suppressed  by  the 
uniqueness  condition  are  marked  accordingly.   To  facilitate 
computation  only  the  restricted  strictly  two  dimensional  problem 
was  considered.   It  suffices  to  consider  the  solutions  to  the 
homogeneous  Integral  equations,  i.e.  functions  representing  flow 
past  a  flat  body,  since  these  solutions  possess  the  same 
singularities  as  do  solutions  of  the  inhomogeneous  problem.   The 
leading  and  trailing  edges  may  be  treated  separately,  so  that 
the  body  becomes  the  positive  or  the  negative  real  axis  with  the 
origin  at  the  leading  or  the  trailing  edge,  respectively.   Under 
these  assumptions  the  flow  becomes  a  similarity  flow:   the  flow 
direction  is  the  same  at  all  points  along  a  ray  from  the  origin. 
A  series  of  flows  were  computed  numerically  for  a  range  of  values 
of  the  parameters  both  in  the   m  >  ■_-  and  m  <  —  regions. 
Figures  16  to  23  show  the  shape  of  the  streamlines  of  these  flows 
near  the  leading  and  trailing  edges. 

The  point  singularity  is  taken  in  the  form 

(^3)  Z2  =  W]_  +  iW2  =  (1  +  iy)  Z~n 

where  n  =  1  -  m  or  m  at  a  leading  or  trailing  edge.   The 
coefficient  y     is  determined  in  terms  of  the  discontinuities 
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across  the  body  in   w\   and  W? .   The  complete  details  of  the 
derivation  of  the  flows  computed  in  this  section  are  given  in 
Appendix  C.   The  form  (43)  taken  from  Z?   may  be  multiplied  by 
a  constant,  thus  affecting  the  flow  variable.   However,  the 
ratio  v/u,  giving  the  streamline  shape,  is  invariant.   The  only 
question  outstanding  is  the  positive  or  negative  direction  of 
the  flow  along  the  streamlines.   The  figures  16  to  23  must 
therefore  be  regarded  as  one  of  two  possible  flows,  the  second 
being  the  reversal  of  the  one  shown. 

The  flow  quantities  have  infinite  values  along  the  char- 
acteristics from  the  leading  and  trailing  edges,  but  the  ratio 
v/u   remains  finite  as  these  characteristics  are  approached. 
However,  there  is  a  discontinuous  change  in  the  streamline  shape 
of  the  characteristic.   In  fact,  it  can  be  seen  from  the  dia- 
grams that  these  singular  characteristics  in  a  sense  may  act  as 
line  sources  or  sinks.   That  is,  although  the  solutions  do 
solve  the  linearized  equations  in  a  generalized  sense,  the 
linearized  equations  are  a  bad  approximation  to  the  non-linear 
ones  in  the  neighborhood  of  the  characteristics.   Thus,  the 
linearized  equations  represent  there  a  bad  approximation  to  mass 
conservation.   It  is  possible  that  the  introduction  of  shocks 
emanating  from  the  leading  and  trailing  edges  may  be  needed. 
The  regional  effect  of  the  line  singularities  of  the  linearized 
solution  may  provide  a  first  approximation  to  the  shock  solution. 


-  33  - 


The  flow  diagrams  may  be  interpreted  as  the  hyperbolic  flow 
carried  along  the  characteristic  superposed  on  the  flow  in  the 
whole  region  near  an  edge  due  to  the  elliptic  point  singularity. 
The  flow  in  the  domain  of  dependence  near  an  edge  characteristic 
is  dominated  by  the  hyperbolic  flow,  but  this  influence  weakens 
away  from  these  characteristics.   At  points  where  the  two 
constituent  flows  are  opposing,  asymptotes  develop.   The  effect 
of  the  point  singularity  outside  the  domain  of  dependence  is 
evident  as  a  mixture  of  a  circulatory  flow  around  an  edge  and  a 
flow  developing  asymptotes. 
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Appendix  A. 

The  equations  leading  to  the  derivation  of  the  integral 
equation  in  the  hyperliptic  restricted  problem  are  set  out. 
The  real  Riemann  invariant  is  given  by  (l6),  which,  with  b  =  0, 
may  be  written  as   r-P  where 


(Al)       P  =  (-|,  v,  £,  t\), 

m 
o 

(  ? 

(A2)      r(A)  =  >  (1-m  )cos  9   +   A  sin  9,      A  cos  9    -    sin 

2         2 
m         m       i 
o      -v  /   o 


i  -  — ^  >  M— ^  ~  i) 

1+A^       1+A^ 


A    is  taken  to  correspond  to  the  characteristic  defining  the 
domain  of  dependence  in  the  lower  half-plane  and   A^   in  the 
upper.   Let   c   be   r(A  )   and   d   be   r(A^).   The  complex 
Riemann  invariant   Zp   is  obtained  from  (Al)  with   A  =  s + it , 
and  may  be  expressed  in  terms  of  its  real  and  imaginary  parts  as 


(A3)      Z2  -  W1  +  W2  =  (a  +  ib)-P 

where 

f     2 
(A4a)     a  =  <  (l-m  )cos  9   +  s  sin  9,      s  cos  9    -    sin  9, 

l-m2(l  +  s2  -  t2)/D,    s[m2(l+  s2+  t2  )  -1  ]/D 


(A4b)      b  =  ft  sin  9,      t  cos  9,      2st  m2D, 
t[m2(l-  s2-  t2)-l]/D  ] 
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(A4c)  D  =  (1  +  s2  -  t2)  +  ds2t2  . 

By  the  elimination  procedures  described  one  integral  equation 
for  the  unknown  function  P(x)   may  "be  obtained: 

(A5)      M-  P(x)  +^  /  f^1   dp  =  k(x) 

where 

(A6a)     [i.   =  (a1b-5-a^b1)(c1di|+ci|d1)  +  (a^-a^  )  (  c^+c^  ) 

+  2c1d1(a^bi|-ai|b    ) 

(A6b)  v   =  (c3d1-c1d3)(a1a4+b1b4)    -    (c^-c^  ^a^+b^  ) 

-  (a2   +  b2)(c3d4-c1+d;5) 

and     k(x)      is  a  known   function. 
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Appendix  B 

Details  of  the  completed  integral  equation  in  the  general 
two-dimensional  problem  are  given.   The  coefficients  of  the 
invariants  associated  with  the  quadratic  factor  are  renamed  to 
avoid  confusion: 


(Bla)  r+  =  e 


(Bib)  r_  =  f 

The  coefficients  of  the  other  two  real  invariants  are  also 
renamed: 


(B2a)  r,  =  cos  j^c 


(B2b)  r^  =  cos  j^d 


and  for  the  complex  invariants 


^  c-oTX,  =  ^  +  ^> 


where   a   and  b   are  real.   It  is  easily  verified  that   a   and 
b   are  given  by  : 
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(B4a)     a  =   M2  s(l+s2+t2)/D  -  sin  0  (cos  0  +  sin  0), 


l-IVr(l+s  -tr)/D,  [(l+s  -tr)(s  cos  0  -  sin  0) 


+  2st2cos  0]/D,  -  sa   +  tb  ,  bs(l+s2+t2 )/D,    b  sin 


(B4b)     b  =  j   M2  t  (l-s2-t2)/D  -  t  sin20,  2st/D, 


[t  ccs  0(l-s2-t2)  +  2st  sin  0]/D,  -  (ta,+  s  b  ), 


bt(l-s2-t2)/D,  0  I 
The  matrices  A,  B  of  (22)  are  given  by 


All  =  "2Mic6dl  +  M2a6bl  -  M5a6dl  +  2M4e6fl 

A12  =  2M1d]_(a1c6-a6c1)/b1-  M^^  -  2M4f -^a^-a^ )/b;L 

A21  =  -2N1c6d1  +  M2a6fl  -  N.a^  +   SN^ 

(B5)  A22  =  AN1d1(a1c6-a6c1)/b1-  2N^f1(a1e6-  a^)/^ 

Bll  =  -  M2aia6 


B 


12  =  a6(M3dl  -  M2V 


B21  =  ° 


B22  =  a6<Vl  "  M2V 
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The  coefficients   M.   in  the  above  may  be  obtained  from  determi- 
nants of  the  matrix 

c3di-  cxd3    c4dr  cxd4    a6(dl-  cx) 


(b6)    m 


a^bn  -  a,b^,    a,,bn  -  anbh     a,-b, 
31    13     41    14     o  1 


a3dl~  ald3    a4di"  aid4    a6'^dl~  al^ 


e3fl"  elf3    e4fl"  elf4   e6fl"  elf6 


V 


after  the  omission  of  the  row  indicated  by  the  subscript.   The 
matrix  N   is  obtained  from  M  by  the  replacement  of  the  second 
row  of   M  by 


Va^f>1"  aif3  '   a4 


fl  *  aif4  '   a6fl"  alf6) 
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Appendix  C. 

The  flows  due  to  the  point  and  line  singularities  near  the 
leading  and  trailing  edges  are  evaluated   in  the  restricted 
problem.   With  the  notation  of  Appendix  A,  the  flow  variables 
are  obtained  from  the  Riemann  invariants  from  the  equation 

4 

(Cl)  P  =  R_1Z  =  /__   F.  .  Z  .   , 

J=l  ±J       J 

where 


(C2)  Z  -  (Z1,  W1,  W2,  Z4) 


and  where   R   is  the   4  X  4  matrix  whose  rows  are  the  vectors 
c_,    a,    b,  _d. 

With  the  homogeneous  boundary  conditions   v,.  =  v,  =  0  and 


u 


b 


the  usual  conditions  on  the  real  invariants  it  follows  that 

(05a)   (z4)t  -  -     *WVt +  Wt    /f24 

(C3b)     (Zl)b  =  -   f22(Wl)b  +  f2J(W2)b   /f2„ 

The  real  invariants  in  the  flow  are  given  by  their  value  traced 
back  along  a  characteristic  to  the  body.   The  flows  for  the 
leading  and  trailing  edges  are  calculated  separately  and  it  is 
convenient  to  take  each  edge  in  turn  as  the  origin  of  co- 
ordinates.  Hence,  the  invariants  have  the  associated  values  at 
any  point  in  the  flow  given  by 
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(C4a)     Z1(x,y)  =H(-y)  H(  -x-T^y)  i  Z±  (x+^y) 


(C4b) 


A 


Z4(x,y)  -  H(y)  H(-x-A^y)  jz4(x+A4y) 


in  the 

trailing  edge 
calculation 


(C5a) 


Z1(x,y)  -  H(-y)  H(x+A4y)  jz^x+A^y) 


(C5b)      Z4(x,y)  =  H(y)  H(x+A4y)   Z4(x+A4y) 


>  t 


in  the 

leading  edge 
calculation 


where  the  body  values  of  the  invariants  on  the  right-hand  side 
are  obtained  from  (CJ>)    and  H(x)  denotes  the  Heaviside  function, 

The  point  singularity  in  the  complex  Riemann  invariant 
according  to  (4j3)  is 


-n 


(C6a)     W-,  =  r   (cos  n^  +  A  sin  nVO 


-n , 


(C6b)      Wp  =  r  "(A  cos  n't    -  sin  nil 


where   r  and  ii     are  defined  by 


(C7) 


z=X+iY=x+sy+    ity   =   re 


111 


The  constant   A   is  determined  from  the  discontinuities  in  W- 
and  Wp   across  the  body.   From  the  relation  between  F(x)  and 
G(x)  and  the  homogeneous  boundary  condition  on   v   it  follows 
that 
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(C8)      b1<(W1)t  -  (Wl)b]  =  a   (W2)t  -  (W2)b 


Hence,  from  (C7),  for  a  trailing  edge  calculation 
(Vt  =  tt,    i^   =    -it) 

(C9)      Y  =  -a1/b1 

and  for  a  leading  edge  calculation  (t/<   =  0,  1>,     =  2ir) 

-a  sin  27m  +  b  (l-cos  27m) 

(CIO)      y  =  -^ - . 

a-,  (l-cos  2-jrn)  +  b,sin  2-rm 

The  flow  variable  vector  P   may  now  be  obtained  from  the 
above  equations.   A  point 

x  =  k(^)cos  j6        y  =  k(j6)sln  j6 
is  chosen,  and  thus 

i/   =  arctan^t  sin  ^/(cos  j6  +   s  sin  jb)\ 

r 


k  (cos  j6  +   s  sin  j6)      +   t  sin' 


The  real  and  imaginary  invariants  are  evaluated  and  substituted 
in  (CI)  for  P.   The  variables   £   and  v   are  given  directly 
and  u   is  obtained  from  the  relation  (3).   The  ratio  v/u  for 
the  streamline  diagrams  discussed  in  Section  7  is  seen  to  be 
independent  of   k  and   r  and  is  a  function  of   6   only.   A 
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similarity  flow  is  thus  obtained 


Appendix  D, 

Specific  solutions  for  flow  quantities  in  the  6   =   „-  case 
are  given  in  this  appendix.   From  the  symmetry  properties  for 
0   =  p-  of  the  differential  equation  and  boundary  condition  it  is 
easily  seen  that  if  the  body  is  an ti symmetric  (i.e.,  (v),  =  (v),), 
and  if  (a,v,w,£,  ,r\  ,C )  (x.y )  is  a  solution  of  the  boundary  value 
problem,  then  (  -a  ,v  ,w,£, ,  -t\  ,  -£ )  (x,  -y )   is  also  a  solution.   Hence 
v,  w,  and  I  are  even  functions  of  y  and   a,  r\,    and  £  are  odd 
in  y.   Similarly,  if  (v),  =  (-v),,   v,  w,  and  £,  are  odd  in  y 
and   c,  i],  £  are  even  in   y.   It  is  evident  that  in  general, 
the  solution  can  be  represented  as  the  sum  of  two  solutions, 
one  depending  only  on  the  even  part,  —  [(v),  +  (v),],  of  the 
assigned  values  of  v,  and  the  other  only  on  the  odd  part, 
o"  [  (v+. )  ~  (v)hJ.   This  decomposition  corresponds  to  the  diagon- 
alization  of  the  system  (22). 

An  explicit  representation  for  boundary  values  of  the 
solution  is  given  explicitly.   It  is  supposed  that  the  boundary 
values  of  v   are  decomposed  as  a  sum 


(Dl)       (v)t  =  g  +  k 


(D2)       (v)b  =  g  -  k  , 
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the  solution  is  a  sum  of  two  functions.   One,  depending  on   g 
alone,  represents  the  flow  about  a  body  with  no  thickness 
(cambered  body)  the  other,  depending  on  R  alone,  represents  the 
flow  about  a  symmetric  body. 

In  the  first  case  (cambered  body)  solution  for  boundary  values 
of  the  density   a  is  given  by 


(D5)      (o)t   =   -(a) 


II! 


.2 

o  (  /     \      v  /  -,  ,  \m-l , .   \-m 
M-S(x)-  -  (1+x)    (1-x) 


r 


(l+s)1"m(l-s)m(s-x)-1g(s)ds 


2   2  ^f*e.\-»-y   u 

/  -,  ,  \m-l  /  ,   \—m 
+  c(l+x)    (1-x) 

where 

(D4)      m  =  m2  =  1  +  £  tan"1  J        |  <  m2  <  1 

(D5a)   u.  =  7s1(l+A2-M2)(l+b2-b2t2AA^+t2)-1(l+b2+b2M2)-1 

+  M  b2(l+b2+b2M2)_1 
o   v        o 

(D5b)      v  =  -t(l+b2+b2A2)(t2+M2-l)(A2+t2)"1(l+b2+b2M2)-1 

and   A,  ,   -A..  ,   it,   -  it  are  the  characteristic  roots  of  the 
linearized  flow  equations. 
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Since   m  >  1/2  ,  the  uniqueness  condition  is  applied  at  the 
trailing  edge.   The  pressure  difference  across  the  body  is  given 

by 


:C6)  (p)   -  (p)   =  2p  u 


2 


(a). 

o 

The  solution  (D3)  in  the  case   b  =  0  with  a  Kutta  condition 
applied  at  the  trailing  edge  corresponds  to  the  solution  obtained 
by  McCune  using  a  Prandtl-Glauert  transformation  in  conjuction 
with  the  incompressible  theory  of  Sears. 

The  singularity  in  the  solution  (D3)  tends  to  the  ordinary 
fluid  dynamics  singularity   (m  =  1/2 )   as  the  magnetic  field 
tends  to  zero. 

As  an  example,  it  is  instructive  to  look  at  the  solution  (D3) 
for  the  flow  past  a  flat  plate  at  the  incidence  a  ,  applying  the 
uniqueness  condition  at  the  trailing  edge.   The  solution  gives 


(D7)     (pt  -  Pb)  =  2pouQa(-v)(\i   +v  )   (sin  irm2 )   f  j 


l-mQ 
1-x     d 


-x, 


The  pressure  singularity  near  the  leading  edge  is  reduced  by  the 
magnetic  field  since   m  >  —  . 

For  the  symmetric  body,  the  solution  for  boundary  values  of 
the  fluid  component   w   of  velocity  is  given  by 
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D(8)  (w),     -    -    (w) 


.1 

-=^  L(X)+^    (l+x)m-1(l+x)-m|(l+s)1-m(l-s)m(s-x)-1K(s)ds, 

d.  .      c.  X.  II 

V1+M-1    _  _J 


m-1, 


+  c1(l+x)m-J-(l-x) 


■m 


where 


(D9a) 


ll,    =    A.M-1    b"1    +  b(M2+t2)(l+^)(t2+A?)-1 
Kl  lo  vo         /v         1 y  1 


(D9b) 
(D9c) 


»!    -    bA1t-1(A2-M2)(l-t2)(t2+A2)-1       , 


1  -1    Vl 

m  =  m,    =  —  tan   "■"  — 

1        ir  (j<1 


and      0  <   m,    <  -= 


(D9d) 


K(x)    =    A,b    1   M~2(l+b2+b2M2)C 
v      '  1  O  O       i 


+   /(M2+t2)(l+A2-M2)(t2+A2)-1-    A^"1         k(x) 


+  ^    (A2-M2)(l-M2-t2)(t2+A2)-1|k(0(^-x)-1    &i, 


-1 


and   C    is  the  constant  value  of  t,      on  the  body.   Here  the 
Kutta  condition  is  applied  at  the  leading  edge. 
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Remark 

The  two  previous  examples  of  the  cambered  and  thickened 
bodies  illustrate  the  general  case  in  which  0  <  m  <  m  <  1  , 
and  one  singularity  is  left  at  each  edge  of  the  general  body 
on  application  of  the  uniqueness  condition.   The  point  on  the 
perpendicular  axis  to   B   at  which  the  singularities  become 
equal  is  given  by 

(D10)  u^v  +  |j,v   =  0 
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Figure  I:    SECTION  OF  THE  RAY  CONE  FOR  A  TWO 

DIMENSIONAL  MAGNETOHYDRODYNAMIC   WAVE. 
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Figure  2    EXPONENTS  m  OF  SINGULARITIES   FOR 
Q  =  IO. 
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Figure  3'-  EXPONENTS  m  OF  SINGULARITIES   FOR    a  =  2 
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Figure  4:  EXPONENTS  m  OF  SINGULARITIES  FOR    OE  =  I.  II 
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a  *  -f  =  0.9 


Figure  5-  EXPONENTS  m  OF  SINGULARITIES    FOR 
a  =  0.9. 
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«  --f  =  0.5 

A? 


Figure  6:  EXPONENTS  m  OF  SINGULARITIES  FOR  a  =  0.5. 
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Figure  7:  EXPONENTS  m  OF  SINGULARITIES 
FOR    a=  O.t. 
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Figure  8:    BEHAVIOR  OF  EXPONENENTS  FOR  SMALL  b. 
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Figure  9:   CRITICAL  FOR   a  =  0.5      b  =  0.5 
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Figure  10-   CRITICAL   CURVES  FOR   a=0.5     b=IO 
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Figure  II     CRITICAL  CURVES  FOR  (2=0.5     b=  1.25 
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Figure  12-    CRITICAL   CURVES   FOR  a  =  0.5    b  =  IO.O 
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Figure  13:    CRITICAL  CURVES  FOR    d  =  2.0    b=  0.5 


-  6l   - 


a  =  2.o 
b»2.o 


1.0  .90 


0.3 


Figure    14:   CRITICAL    CURVES  FOR    a  =  2.0    b=2.0 
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a  =  2.0 
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Figure  15:    CRITICAL  CURVES  FOR   a  =  2.0     b  =  10. 0 
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Figure  16:      sin#=l         MQ  =  I  a  =  .5 

LEADING   EDGE  FLOW      |-m 


.022584 
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Figure  17.      sin0=l  M-l.  a  =  .5 

TRAILING  EDGE  FLOW         m=  .977416 
FLOW  SUPPRESSED 
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Figure  18 


sin  a  =  .15  M0  =  .45 

LEADING   EDGE  FLOW 
FLOW  SUPPRESSED 


a  =  2 

I  -m  = 


68842 
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Figure  19: 


TRAILING  EDGE  FLOW 
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Figure  20:    sin0=.l5     M0=2.2 

LEADING  EDGE  FLOW 
FLOW  SUPPRESSED 


a  =  10. 

l-m=  .53228 
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Figure  21:     sin#=.l5  Mq=  2.2        a-    10. 

TRAILING  EDGE  FLOW       m  =  .46772 
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Figure  22:  sin#  =  .45  M0=  .95 

LEADING  EDGE  FLOW 
FLOW   SUPPRESSED 


CL  =  .5 
l-m  =  .62326 
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Figure  23:    sin0=45  M0=  .95  a  =.5 

TRAILING  EDGE   FLOW       m  =  . 37674 
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